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Abstract 

By using first-principles calculation, we have found that a family of 2D transition metal dichalco¬ 
genide haeckelites with square-octagonal lattice MX2-4-8 (M—Mo, W and X=S, Se and Te) can 
host quantum spin hall effect. The phonon spectra indicate that they are dynamically stable and 
the largest band gap is predicted to be around 54 meV, higher than room temperature. These will 
pave the way to potential applications of topological insulators. We have also established a sim¬ 
ple tight-binding model on a square-like lattice to achieve topological nontrivial quantum states, 
which extends the study from honeycomb lattice to square-like lattice and broads the potential 
topological material system greatly. 
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I. INTRODUCTION 


Inspired by the impressed progress in theory and applications, nnmerous researchers tnrn 
their attentions to two-dimensional (2D) systems especially miraculous graphen^. Recently, 
another 2D system, transition metal dichalcogenides (TMDs) MX 2 with M=Mo, W, Ti, etc. 
and X=S, Se, Te, has been widely explored due to wide range of electronic properties!^ and 
easy fabrication. Bulk TMDs are composed of 2D X-M-X layers stacked on top of each other. 
The bonding within those trilayer sheets is covalent while the coupling between adjacent 
sheets is weak van der Waals (vdW) interaction. 2D TMD layers can be manufactured not 
only by mechanicaP^ and chemical exfoliatiorpESl of their layered bulk counterparts, but 
also by chemical vapor deposition (CVD)P! or two-step thermolysiJ^. However, defect will 
be inevitably produced during the manufacture, and it will modify the electronic structure 
significantly in low dimensional system. 

Graphene with 5-7 defects, which is often called Haeckelites in honor of the German 
biologist and naturalist Ernst Haeckel, has been theoretically proposed about twenty years 
agcJEEl. In contrast to comprehensive understanding of the defect in graphen^^^, defect 
in 2D TMDs may just launch on. There are still so many issues need to be tackled. W. 
Li et al^ and H. Terrones et alW^ proposed that a new planar sheet could be generated 
from original hexagonal TMDs when introducing 4-8 defects. This new planar sheet is 
TMD Haeckelites with square-octagonal lattice. The periodic 4-8 defects have been ob¬ 
served in the grain boundaries of M 0 S 2 and most probably will exist in other TMD^, too. 
High-resolution transmission electron microscopy (HR-TEM) is a powerful instrument to 
selectively suppress or enhance bond rotations and produce defects in sample due to the 
ballistic procedure between high energetic electrons and sample atoms. Using aberration 
corrected HR-TEM device, a disordered graphene Haeckelite has been produced in sitrP^. 
Combing first-principles calculations and HR-TEM experiments, H. Komsa demonstrated 
that it is possible to observe defect formation under exposure to an 80 keV electron beam 
in M 0 S 2 systerrpSI. 

Recently, Qian et al.!^ predicted that some 2D TMDs with IT' structure can be large-gap 
2D topological insulators (TIs) though most of them are in IH structure and are not TIs. It is 
natural to ask whether TI state can exist in the defected 2D TMDs based on IH structure. 2D 
TP^ES^ also known as quantum spin hall (QSH) insulator, was firstly proposed in graphene. 
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where spin-orbit coupling (SOC) opens a band gap at the Dirac point. On account of weak 
SOC strength, the band gap is so small (order of 10“^ meV)p2 that this proposal is hardly 
to be verified by experiments. So far, QSH effect has only been observed in HgTe/CdT^ 
and InAs/GaStP^ quantum wells. Both of them require precisely controlled MBE growth 
and ultralow temperature. The study of 2D TI has been seriously hampered due to lack of 
proper materials with large band gap, stable structure and easy fabrication.!^ In this work, 
based on first-principles calculations, we find monolayer of WX 2 and M 0 X 2 Haeckelite are 
2D TIs and the largest band gap is around 54 meV. Distinguished from other predicted 
QSH materiald^^!® based on honeycomb lattice, these MX2-4-8 Haeckelites have square- 
like lattice and a simple tight-binding model with one orbital per site and four sites per 
unit cell has been established to achieve topologically nontrivial QSH state. Such extension 
from honeycomb lattice to square-like lattice have largely broad the potential candidates for 
topological materialsl ^^ * ^^ * 

The paper is arranged as follows. In section H we will introduce the details of first- 
principles calculations. In section HI, the calculation results are presented and TB analysis 
is performed. Finally, section IV contains a conclusion of this work. 


II. CALCULATION METHOD AND CRYSTAL STRUCTURE 

First-principles calculations were carried out using the projector augmented wave methocpSEl 
implemented in Vienna ab initio simulation package (VASP)^^^*^. Exchange and correlation 
potential was treated within the generalized gradient approximation (GGA) of Perdew- 
Burke-Ernzerhof typ^^. SOG was taken into account by the second variation method 
self-consistently. The cut-off energy for plane wave expansion was 500 eV. The k-points 
sampling grid in the self-consistent process was 9x9x3. The crystal structures have 
been fully relaxed until the residual forces on each atom were less than 0.001 eV/A. The 
crystal parameters for all TMD Haeckelites are shown in Table |i| A vacuum of 20 A between 
layers was considered in order to minimize image interactions from the periodic boundary 
condition. PHONOPY has been employed to calculate the phonon dispersiorP^. To explore 
the edge states of TMD Haeckelites, maximally localized wannier functions (MLWFs) for 
the d orbitals of W and p orbitals of S have been constructed and used to get ab initio tight- 
binding (TB) hamiltoniarP^. Atomic SOC is added to the TB hamiltonian by fitting the 
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FIG. 1: (color online). Top view(a) and side view (b) of relaxed WS2-4-8. Gray ball is W 
and yellow ball is S. The primitive cell is shown in light blue rectangle, (c) 2D and 
projected edge first BZ with high symmetry point (red dots), (d) Phonon dispersion of 

WS 2 - 4 - 8 . 


first-principle calcualtions. The projected edge states were obtained from the TB through 
an iterative methocP^®®. 

All TMD Haeckelites have the same non-symmorphic space group Pbam (D®/^). Except 
the difference of lattice constants, all other TMD Haeckelites have nearly the same properties 
as WS 2 . So we choose WS 2 - 4-8 later as an example and the results of other TMD Haeckelites 
can be found in the Appendix. The relaxed crystal structure and Brillouin zone (BZ) for 
WS 2 - 4-8 are shown in Fig. It is noted that bond length between W 1 and 3 is reduced 
compared to that in honeycomb lattice. As we will see later, this bond is vital to topological 
phase transition. The dynamic stability of WS 2 - 4-8 has been investigated by calculating 
it phonon spectrum. Imaginary frequencies can not be found in the phonon dispersion of 
WS 2 - 4 - 8 , which indicates the structure is stable (Fig. [^d)). 
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TABLE I: Lattice parameters and band gaps for some TMD Haeckelites MX2-4-8. 


Structure 

a(A) 

b(A) 

Gap (meV) 

WS 2 

6.34 

6.41 

53.82 

WSe2 

6.40 

6.86 

30.03 

WTes 

6.65 

7.38 

14.743 

M 0 S 2 

6.36 

6.33 

13.38 

MoSe2 

6.65 

6.59 

26.80 



FIG. 2: (color online). GGA (a) and GGA+SOG (b) band structures of WS2-4-8. The 
calculated edge states for X (c) and Y (d) edge, respectively. 

III. RESULTS 

Both GGA and GGA+SOG band structures of WS2-4-8 are shown in Fig. 2. Along 
F-X direction, there is a band crossing in GGA band structure. The little group of k point 
located on this direction is C 2 V The two crossing bands belong to different irreducible 
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representations and such band crossing is protected by {C' 2 x||| 0 } operation. When SOC 
is included, as we can see in Figj^ (b), WS2-4-8 is a well defined insulator with indirect 
band gap around 53.8 meV. Gaps for other TMD Haeckelites are listed in the Table |Ij 
For system possesses both time reversal (TR) and space inversion symmetry, the parity 
criterion proposed by Fu and Kane is a convenient method to judge its band topology 
Z 2 number.l^ Since the space group of it is non-symmorphic, similar as that in single layer 
ZrTes and HfTe^, all the bands at the three time reversal invariant momenta X, S, Y having 
degeneracy of even and odd states. Only the band inversion at F can leads to nontrivial 
band topology. It is true that the total parity of occupied states at F is -1 and the other 
three are +1. Therefore, the topological index Z 2 equals to 1, which means WS2-4-8 is a 
QSH insulator. Considering the possible underestimation of band gap of GGA, non-local 
Heyd-Scuseria-Ernzerhof (HSE06) hybrid functionaP^ is further supplemented to check the 
topological property. The band topology does not change. 

Since the existence of nontrivial edge states is the hallmark of QSH effect, we have 
calculated the edge states of WS2-4-8. As shown in Fig. (c, d), there is edge Dirac cone 
dispersion connecting the bulk occupied and unoccupied states with Dirac point at F for 
both X and Y edges. 

In order to understand the band inversion process at F point explicitly, a Slater-Koster 
TB has been constructed. Obviously, the main physics comes from the isolated group of four 
bands around the Fermi level. The band character and projected density of states (PDOS) 
analyses indicate that the low energy bands near the fermi level are mainly contributed by 
the and orbitals of W. The can mix with due to the distortion of 4-8 

defects. For simplicity, only ^^2 is taken into account for each W and one primitive cell 
contains four W atoms. Therefore, it is possible and reasonable to construct a 4 x 4 TB 
hamiltonian for non-SOC case. The non-SOC TB hamiltonian with four localized d^^ bases 
\uji) (i=l, 2, 3, 4) can be written as follows in momentum space 


Ho{k) 


^ 0 2ticos{\kx) 

2ticos(^kx) 0 

2tsCos{^ky) 

y 2^3 cos(|fcj^) 


^262^ 

2^3 cos{\ky) 

0 

2ti cos (4^3,) 


2^3 cos{lky) ^ 

2ti cos(4A:3.) 

0 / 


( 1 ) 
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where the hopping parameters are defined as 


ti = {wi\Ho\w2) 

(2) 

t2 = {wi\Ho\w3) 

(3) 

ts = {Wi\Ho\w4). 

(4) 


\ui) means the orbital (1^2 on z-th W atom labelled in Fig. (a). However, one should 
note that periodic boundary condition is sacrificed here in order to get a concise form of Eq. 
Q. When the calculated properties concern global phase such as berry phase, we need to 
transform H^ik) into another form which satisfies periodic boundary condition. 

Similar to graphene, SOC is a second order effect for WS 2 - 4 - 8 . The intrinsic atomic SOC 
for W is of the order of 200 meV while it is about 4 meV for graphene. Compare to the 
extremely small gap for graphene, a large band gap (53.8 meV) is obtained for WS 2 - 4-8 
at last. The hybridization between jo;* t) S’lid \u}j 4-) is zero due to symmetry. If TR 
symmetry is preserved, we will have Hl'lik) = H^{kY'. Therefore, even the hamiltonian 
size will be doubled when SOC is taken into account, we can still focus on spin up (spin 
down) subspace only. Spin down (spin up) subspace can be obtained using above restricted 
conditions. Considering all the symmetries, we obtain a generic matrix form for H^{k). 


hVM = 


^ 0 2Ai cos( 4 A:a,) 0 2A3 cos( 4 A:^) ^ 

2 A 4 cos(4A:a,) 0 2A3Cos(4A:^) 0 

0 2A3Cos(4fcj^) 0 2Ai cos(|A:a;) 

y 2 A 3 cos(4fcj^) 0 2Xlcos(^kx) 0 J 

where Ai and A 3 (pure imaginary numbers) are defined as 


Ai = {wi t \Hs 0 \w 2 t) 
A 3 = {wi t \Hs 0 \W 4 t) 


(5) 


( 6 ) 

(7) 


Only the total sign(tit 2 t 3 ) makes sense. Sign(tit 2 t 3 ) < 0 can be inferred by inverting the 
whole bands of sign(tit 2 t 3 ) > 0 . So we take > 0 , ^2 > 0 and > 0 . As discussed above, 
only the band inversion at T point can change the band topology. The four eigen states and 
their parities can be obtained explicitly as following: 
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FIG. 3: (color online). Non-SOC (black lines) and SOC (red lines) band structure 
evolution calculated with different t 2 in TB model. The parameters are Ai = —A 3 = 0.08 i, 
ti = 1 , G = 0 . 6 , ^2 = 2.2 (a), ^2 = 2.0 (b), ^2 = 1-8 (c), ^2 = 1-2 (d) and ^2 = 1-0 (e). (f) 
ti = 0.6, G = 1 and ^2 = T 8 . Parity information at P has been given. 


El 

= 2 tl + t 2 + 2 t 3 

parity 

+ 

( 8 ) 

E 2 

= - 2 tl + t 2 - 2 t 3 

parity 

+ 

(9) 

E 3 

= 2tl -t2- 2t3 

parity 

— 

( 10 ) 

El 

= -2tl -t2 + 2t3 

parity 

— 

( 11 ) 


Obviously, the Ei has the highest energy. The band inversion happens between E 2 and 
E 3 or E 2 and E 4 can lead to QSH state. That means the band inversion will exist as long 
as |t 2 |<niax(| 2 ti|, | 2 t 3 |). Fig. show the band structure calculated with the TB model 
with different sets of parameters. When |t 2 | > 2|ti| > 2 |t 3 | (Fig. (a)), the bonding and 
anti-bonding states are far away from each other. This is a trivial insulator. There will be a 



















band touching at F point if \t 2 \ = 2|ti| (Fig. |^(b)). When 2|ti| > > 2 1 ^ 31 , band inversion 

will happen (Fig. (c)) and the system enters into QSH state. The two valence bands {E 2 
and E 4 ) with opposite parity will be degenerate at F if |t 2 | = 2|t3| (Fig. (d)). When 
2|ti| > 21 ^ 3 ! > 1 ^ 2 !, they will separate each other (Fig. [^(e)). Topological phase transition 
will not take place in this parametric region from Fig. |^c) to (e), since the band inversion 
process occurs between two occupied valence bands. Another interesting thing is that ti and 
ts control the position of band crossing in non-SOC band structure. If > ts, the crossing 
will be situated at F — X direction (Fig. (c)) and it will move to F — Y direction (Mo case 
in Appendix) for ti < (Fig. i(f)). 

Uniaxial strain effect has also been investigated with first-principles calculations, which 
can be used to tune the relative strength of these hopping parameters and control the 
topological quantum state transition. If we fix crystal constant a and decrease b, will 
increase faster than ti and t 2 - This means the band crossing in non-SOC band structure will 
move from F — X to F — Y and uniaxial strain is always beneficial for band inversion. Similar 
results can be get for uniaxial strain along b-axis. Therefore the topology in WS2-4-8 is 
robust against uniaxial strain. 


IV. CONCLUSION 

In summary, we have performed first-principles calculations for the electronic properties 
of TMD Haeckelites MX2-4-8 and found that they are 2D TIs. A simple TB model with 
one orbital per site and four sites per unit cell has been established to understand the band 
inversion mechanism. The nearest hopping parameter t 2 is vital to trigger the topological 
phase transition and can be tuned through lattice strain effect. Such simple square-like 
lattice model to achieve various topological quantum states is very stimulating. It will lead 
further study based on square lattice instead of honeycomb lattice and largely extend the 
searching range for topological materials. 
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V. APPENDIX: OTHER TMD HAECKELITES MX2-4-8 



FIG. 4: (color online). Top views of crystal strnctnre of WSe2-4-8 (a), WTe2-4-8 (a), 

M 0 S 2 - 4-8 (c) and MoSe2-4-8 (d). 


13 






















FIG. 7: (color online). Edge states for WSe2-4-8 (a, b), WTe2-4-8 (c, d), MoS2-4-8 (e, f) 
and MoSe2-4-8 (g, h). Left panel for X edge and right panel for Y edge. 
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